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Abstract 

Background: Cotton worm is one of the main insects of soybean in southern China. Plants may acquire defense 
mechanisms that confer protection from predation by herbivores. Induced responses can lead to increased 
resistance against herbivores in many species. This study focuses on searching changed proteins in soybean 
defense response induced by cotton worm feeding. 

Results: Ten protein spots that are changed in abundance in response to cotton worm feeding were identified by 
Two-dimensional gel electrophoresis (2-DE). A total of 11 unique proteins from these spots were identified by 
MALDI-TOF MS. The mRNA and protein relative expression levels of most changed proteins were up-regulated. 
These proteins were mainly involved in physiological processes, including active oxygen removal, defense signal 
transduction, and metabolism regulation. 

Conclusion: This is the first proteomic analysis of the soybean defense response induced by cotton worm. The 
differentially expressed proteins could work together to play a major role in the induced defense response. PAL 
and SAMS were up-regulated at both the protein and mRNA levels. These genes can be strongest candidates for 
further functional research. 
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Background 

As an important economic crop, soybean provides sig- 
nificant sources of fatty acids, proteins, vitamins, miner- 
als and other nutrients for humans and animals, and it 
also has nonfood uses, such as in the production of 
industrial feedstocks and combustible fuels [1,2]. How- 
ever, insect pests can adversely affect the yield and qual- 
ity of soybean. Cotton worm is one of the main insect 
pests of soybean in southern China [3,4]. Therefore, 
improvement in resistance to cotton worm is one of the 
main soybean breeding objectives. 

During the process of evolution, plants have acquired 
defense mechanisms that confer protection from preda- 
tion by herbivores. Two modes of resistance to herbi- 
vores exist: constitutive resistance, which is expressed 
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independent of an attack; and induced resistance, which 
is activated only after the plant is attacked or otherwise 
injured [5]. Compared with constitutive resistance, 
induced resistance may be more durable and compro- 
mises plant fitness less by either decreasing further her- 
bivore damage or increasing plant tolerance to herbivory 
[6]. Induced responses that lead to increased resistance 
against herbivores have been reported for over 100 spe- 
cies of plants, such as Ambidopsis [7-9], tobacco [10], 
tomato [11,12], rice [13], and soybean [14-16]. The 
resistance responses induced in plants change with dif- 
ferent attackers. For example, the level of resistance 
induced by soybean looper herbivory with subsequent 
bean leaf beetle feeding was higher than that induced by 
bean leaf beetle herbivory with subsequent soybean 
looper feeding [16]. Resistance induced by the soybean 
looper, Mexican bean beetle and corn earworm in soy- 
bean has been reported in many studies [14-17]. 
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However, induction of resistance in soybean by cotton 
worm has not been reported. 

Many studies have investigated the induced defense 
response mechanism of plants. The plant hormones sal- 
icylic acid (SA), jasmonic acid (JA) and ethylene (ET) 
are the main players in the regulation of signaling net- 
works involved in induced defense [18-21]. In general, 
pathogens are more sensitive to SA-dependent 
responses, whereas herbivorous insects and necrotrophic 
pathogens are resisted by JA/ET-dependent defenses 
[22-24]. Ample evidence indicates JA is the main signal- 
ing molecule that mediates a plant's defense system 
against herbivores [19,23]. A classic example is the 
observation that following attack by Manduca sexta lar- 
vae, tomato {Solanum lycopersicum) leaves accumulate 
JA, which results in the activation of genes encoding 
proteinase inhibitor proteins that inhibit digestive Ser 
proteinases of herbivorous insects and reduce further 
insect feeding [8,25]. Induced defense response mechan- 
isms also involve the production of defensive com- 
pounds, such as proteinase inhibitors that affect insect 
feeding [25,26], volatile organic compounds that attract 
parasitoids and predators of the herbivores that feed on 
the plant [18,27] and extrafloral nectar that arrests car- 
nivorous arthropods on herbivore-infested plants [28]. 

Proteomics is a powerful tool for the study of protein 
dynamics, especially in plant stress responses [29,30]. 
Analysis of protein profiling is important to under- 
stand how genes/proteins are regulated. In the present 
study, we found that soybean resistance to cotton 
worm increased after induction of resistance by cotton 
worm feeding. To investigate which proteins are 
involved in this induced defense response, we used 
proteomic approaches to provide an overview of the 
cotton worm-induced defense response in soybean. 
Ten protein spots showed differential expression pat- 
terns in soybean leaves in response to cotton worm 
feeding, and 11 proteins were identified. Quantitative 
RT-PCR (qRT-PCR) analysis indicated differentially 
expressed proteins also showed changes in mRNA 
expression level. 



Results and discussion 

Induction of soybean defense response by cotton worm 
feeding 

To determine whether cotton worm feeding can induce 
resistance in soybean leaves, we carried out a dual-choice 
test and force-feeding experiment (Figure 1). Compared 
with treated leaves, the cotton worm larvae preferred to 
eat control leaves, where PI < 1 at the four sampling 
timepoints (Figure 1A). The relative growth rate of cot- 
ton worms feeding on control leaves was higher than 
those feeding on treated leaves except at 6 h (Figure IB). 
In the force-feeding experiment, the relative growth rates 
of cotton worms feeding on control leaves and those 
feeding on treated leaves were similar. At this timepoint 
fifth-instar cotton worms were used, which eat greedily 
and were not sensitive to the food source [31]. In the 
dual-choice test and force-feeding experiment, the levels 
of induced resistance were slightly different at the four 
sampling timepoints, but these results still indicated that 
cotton worms can induce resistance in soybean. 

Plants do not exhibit induced resistance to insects unless 
subjected to a certain degree of damage. In pre-experi- 
ments, soybean treated with cotton worm feeding shown 
induced resistance at a damage level of 20-30% and the 
levels of resistance were not increased significantly when 
the damage area increased. The capacity for induced resis- 
tance showed no difference also at a damage level between 
15-45% in mountain birch [32]. However, when leaf 
damage constitutes 50% of the leaf area, induced resistance 
was produced in Quercus garryana [33]. The timing of 
induced resistance to herbivory is known to vary between 
plant species. Usually, in herbaceous plants, the longevity 
of induced resistance is several days [34-36] . The timing of 
induced resistance varied with the insect species feeding 
on soybean. The longevity of resistance induced by C. tri- 
furcata was two weeks [16]. The induced response was 
effective for 3 days after damage against the Mexican bean 
beetle in soybean. After 15 days, the strength of resistance 
had declined and by 20 days, all four genotypes exhibited 
induced susceptibility [37]. In this study, the longevity of 
induced resistance produced by cotton worm was 5 days. 
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Figure 1 Identification of induced resistance. A, Feeding preference indices in dual-choice tests.B, Relative growth rate of cotton worms fed 
with control (CK) and treated FT) leaves. Error bars represent the standard deviation. * Significant (P < 0.05). 
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Figure 2 Representative silver-stained 2D-PAGE gel of separated 
soybean leaf proteins. Proteins were separated in the first dimension 
on a nonlinear IPG strip, pH 3.0-10.0, and in the second dimension on 
a 12% polyacrylamide SDS-gel. Quantitative image analysis revealed a 
total of 10 spots that changed in abundance. 



Changes of proteins in soybean leaves induced by cotton 
worm feeding 

Soybean leaf proteins extracted 0 h, 6 h, 12 h and 24 h after 
cotton worm feeding and from the controls were analyzed 
by 2-DE. Three independent experiments were conducted 
to ensure that the changes in protein abundance at each 
timepoint were reproducible and significant. Software quan- 
tification showed that, although the protein expression pro- 
files between the treated and control samples were similar, 
the abundance of some spots changed significantly with 
cotton worm feeding over time (Figures 2, 3 and 4). Ten 
spots showed changes in abundance in response to cotton 
worm feeding, most of which were up-regulated (Figures 2 
and 3). The different expression patterns of the spots might 
imply different roles for these proteins in plant defense 
responses induced by cotton worm feeding. 

Identification of proteins involved in cotton worm- 
induced defense response in soybean 

All 10 differentially expressed proteins spots were sub- 
jected to in-gel digestion and analyzed by MALDI-TOF- 
MS. The peptide mass fingerprints obtained were used 
to search the NCBI database using Mascot http://www. 
matrixscience.com. There were 11 proteins identified 
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Figure 3 Enlarged maps of the 10 differentially expressed protein spots. CK: control; T: soybean leaves treated with cotton worm feeding. 
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Figure 4 The volume changes of 10 differentially expressed spots. Control (black) and treated (gray). A - Jrepresent spot 1-10 
respectively. Error bars represent the SD. The scales are different for the different proteins. 



from these 10 protein spots (Table 1). Spots 5 and 10 
contained two protein species, respectively. Spots 2 and 
9 were identified as ribulose-1, 5-bisphosphate carboxy- 
lase/oxygenase large subunit. When using 2-DE, co- 
migrating proteins can present a major problem for 
quantitative protein expression comparison [38,39]. In 
this study, real-time RT-PCR was used to help 



determine which proteins showed altered expression. 
Occurrence of the same protein in several spots suppo- 
sedly represents differently modified protein species that 
might be of biological relevance [38]. Interestingly, the 
observed Mr/PI of some proteins was different from 
their theoretical Mr/PI. The reason might result from 
an unusual posttranslational modification. 



Fan ef al. Proteome Science 2012, 10:16 
http://www.proteomesci.eom/content/1 0/1/16 



Page 5 of 1 1 



Table 1 Identification of differentially expressed protein spots by MALDI-TOF MS 



Spot no. 


Protein name 


Mr(KDa)/PI 


Sequence coverage (%) 


Score 


Accession no. 






Theoretical 


Experimental 








1 


serine/threonine kinase-related protein 


14.7/9.4 


42.0/6.0 


-16 


76 


gi|34099884 


2/9 


ribulose-1,5- bisphosphate 


23.3/6.1 


24.0/7.7 


23 


90 


gi| 1 5781 2656 




carboxylase/oxygenase large subunit 


51 .8/6.2 


31 .1/0.81 


30 


1 72 


gi| 14599574 


3 


Vegetative storage protein A 


29 2/8 8 


24 4/6 1 


62 


1 7'\ 


^jl 1 1 T) A 1 A r 

gi 1 3414b 


4 


LyLUSUIIL pi lUSpi IULJ lyLtrl d Lc KllldSfcr 


42.3/5.7 


43.6/6.4 


26 


1 03 


nilQ93n771 
y i y/ju/ / i 


5 


ATP synthase beta subunit 


50.0/5.7 


42.4/6.4 


32 


81 


gi|20269424 


5 


malate dehydrogenase 


36.1/8.2 


42.4/6.4 




68 


gi|5929964 


6 


phenylalanine ammonia-lyase 


46.5/6.0 


40.3/6.7 


56 


347 


gi|81 807 


7 


S-adenosylmethionine synthetase 


42.9/5.4 


47.4/6.1 


25 


88 


gi| 10443981 


8 


predicted protein 


60.3/5.8 


47.0/6.5 


24 


75 


gi| 168040478 


10 


ascorbate peroxidase 


27.1/5.7 


29.8/5.9 


45 


110 


gi| 12229897 


10 


20S proteasome alpha subunit 


27.4/5.8 


29.8/5.9 


41 


104 


gi| 12229897 



We verified the cotton worm-induced defense 
response profile of abundance-changed proteins at 
mRNA level by real-time qRT-PCR. The mRNA levels 
of most proteins were up-regulated in response to cot- 
ton worm feeding (Figures 4 and 5). ATP synthase beta 
subunit and malate dehydrogenase were located in the 
same spot, and the mRNA levels for both proteins were 
up-regulated, which indicated expression of both protein 
species in spot 5 was changed in response to cotton 
worm feeding. Similarly, expression of ascorbate peroxi- 
dase and 20S proteasome alpha subunit in spot 10 both 
changed in response to cotton worm feeding. 

Levels of mRNA are not always consistent with the 
levels of the corresponding proteins [40,41]. Three 
potential reasons for the lack of a strong correlation 
between mRNA and protein expression levels have been 
proposed: (i) translational regulation; Because of many 
complicated and varied post-transcriptional events, the 
transcriptome might not always be reflected at proteome 
levels [42-44]. (ii) differences in protein in vivo half-lives 
[42,43]; and (iii) significant experimental error, there is a 
significant amount of error and noise in both protein 
and mRNA experiments that limit our ability to get a 
clear picture. 

Differentially expressed proteins 

Spot 1 was up-regulated in response to cotton worm 
feeding, which was identified as serine/threonine 
kinase (STK)-related protein (Figure 4A and 5A). As a 
signal protein, STK regulates a variety of cellular pro- 
cesses [45,46]. This study found that STK involved in 
soybean defense responses induced by cotton worm 
feeding. Previously, it has been reported that STK was 
shown to be involved in plant defense responses and 
induced by low temperature, Soybean mosaic virus, 
and SA [47-49]. 



Spots 2 and 9 were up-regulated and were identified 
as ribulose-l,5-bisphosphate carboxylase/oxygenase 
(RubisCO) large subunit(Figures 4B, I and 5B). RubisCO 
is the most abundant protein in plants and plays a 
major role in photosynthesis. RubisCO large subunit is 
up-regulated under stress, which might be indicative of 
the protective and more active action of carbon assimi- 
lation machinery [50]. 

Spot 3 was identified as vegetative storage protein a 
(VSPa). This protein was up-regulated at both protein 
and mRNA levels (Figures 4C and 5C). The best-charac- 
terized VSPs are soybean VSPa and VSPp VSPs are 
thought to serve as a transient reserve that sequesters 
unused amino acids during plant development. Once 
new seed production begins, stored VSPs presumably 
make nitrogen and other nutrients immediately available 
for seed development [51]. This hypothesis, however, 
was not supported by transgenic soybean in which VSP 
expression was abolished and seed production was unaf- 
fected even under nitrogen-deprived conditions [52]. 
Thus VSPs may serve other functions beyond source- 
sink interaction or plant productivity. Arabidopsis VSP 
transcripts are induced by mechanical wounding, JA, 
insect herbivory, and osmotic and nutritional stresses 
[9,19,53-55], which is a common response shared by 
many genes encoding anti-insect proteins. Liu et al. 
found AtVSP2 is an anti-insect acid phosphatase [56]. 
Whether soybean VSPa has an anti-insect function 
needs further study. 

Spot 4 showed up-regulation at 0 h and was identified 
as cytosolic phosphoglycerate kinase (PGK) (Figure 4D). 
It was up-regulated at all sample timepoints at the 
mRNA level (Figure 5D). PGK is an enzyme that cata- 
lyzes the formation of ATP from ADP and vice versa. In 
the second step of the second phase in glycolysis, 1,3- 
bisphosphoglycerate is converted to 3-phosphoglycerate, 
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Figure 5 The expression of 10 differentially expressed spots at the corresponding mRNA level A represents spot I B represents spot 2 
and 9. C - Grepresent spot 3-7 respectively. J represents spot 10. A, serine/threonine kinase-related protein.B, ribulose-1 ,5- bisphosphate 
carboxylase/oxygenase large subunit. C, Vegetative storage protein A. D, cytosolic phosphoglycerate kinase.E- 1, ATP synthase beta subunit. E- 2, 
malate dehydrogenase. F, phenylalanine ammonia-lyase.G, S-adenosylmethionine synthetase..)- 1, ascorbate peroxidase. J- 2, 20S proteasome 
alpha subunit. Error bars represent the SD. The scales are different for the different genes. 



forming one molecule of ATP. If the reverse occurs, one 
molecule of ADP is formed. This reaction is essential in 
most cells for ATP generation in aerobes, for fermenta- 
tion in anaerobes, and for carbon fixation in plants. 
Under stress conditions, do plants need to energy has 
not yet been reported. 



Spot 5 was identified as ATP synthase P subunit and 
malate dehydrogenase (MDH). These proteins were up- 
regulated at the mRNA level (Figure 5E-1 and 5E-2). 
ATP synthase is an important enzyme that generates 
energy for the cell to use through the synthesis of ATP. 
ATP synthase P was up-regulated in the rice defense 
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response induced by probenazole [39]. MDH is an 
enzyme in the citric acid cycle that catalyzes the conver- 
sion of malate into oxaloacetate (using NAD+) and vice 
versa (this is a reversible reaction). MDH is not to be 
confused with malic enzyme, which catalyzes the con- 
version of malate to pyruvate, producing NADPH. 
MDH is involved in many physiological activities, such 
as the organelle shuttle system, plant disease resistance, 
and metabolism of reactive oxygen species. 

Spot 6 was identified as phenylalanine ammonialyase 
(PAL). This protein was up-regulated at both the pro- 
tein and mRNA levels (Figures 4F and 5F). It is clear 
that PAL plays a key role in induced disease resistance 
based on induction of its expression by pathogens or eli- 
citors, and on studies using transgenic approaches in 
tobacco [57,58]. PAL is the first enzyme in the phenyl - 
propanoid pathway. In plants, PAL can supply the pre- 
cursors for lignin, flavonoid pigments, protectants, and 
furanocoumarin phytoalexins [59]. Lignin has a positive 
effect on plant resistance to insects. Based on its up-reg- 
ulation in the soybean defense response induced by cot- 
ton worm feeding, PAL might play a key role in cotton 
worm-induced resistance. 

Spot 7 was identified as S-adenosylmethionine synthe- 
tase (SAMS). This protein was up-regulated at both the 
protein and mRNA levels (Figure 4G and 5G). Previous 
studies showed that increased expression of SAMS was 
induced by mechanical wounding, JA, insect herbivory, 
chilling, and salinity stress [60-63]. SAMS is a key 
enzyme that catalyzes the formation of S-adenosyl- 
methionine (SAM) from methionine and ATP. SAM is a 
precursor of polyamines and ET. Polyamines play 
important roles in plant defense to a variety of environ- 
mental stresses [64,65]. ET as a plant hormone that 
plays an important role in induced defense. Thus SAMS 
might play a key role in soybean induced resistance and 
warrants further study. 

Spot 10 was identified as ascorbate peroxidase (APX) 
and 20S proteasome a subunit. These proteins were up- 
regulated at the mRNA level (Figure 5J-1 and 5J-2). APX 
is an enzyme in the ascorbate-glutathione cycle. The 
ascorbate-glutathione cycle is an important antioxidant 
protection system against H 2 0 2 generated in different cell 
compartments. In this study, APX may contribute to 
induced resistance. Overexpression of the APX gene not 
only enhances the plant's antioxidant capacity, but also 
can enhance plant tolerance to drought, heat, salt and 
pathogens [66,67]. In eukaryotic cells, most proteins in the 
cytosol and nucleus are degraded via the ubiquitin-protea- 
some pathway. The 26S proteasome is responsible for the 
degradation of damaged and misfolded proteins. A pre- 
vious study showed that elicitation of defense reactions in 
tobacco cells by cryptogein, a proteinaceous elicitor of 



plant defense reactions, leads to a rapid and differential 
accumulation of transcripts corresponding to genes that 
encode defense-induced subunits of 20S proteasome [68] . 

In summary, all of the above proteins showed changed 
abundances induced by cotton worm feeding on soybean 
leaves. Thus the induced resistance might be a result of 
the combined effect of these proteins. 

Conclusions 

In order to better understand the defense response 
induced in soybean, we investigated the changes in the 
soluble proteome of soybean leaves damaged by cotton 
worm feeding using 2-DE and MALDI-TOF-MS. A total 
of 10 spots changed in abundance, from which 11 pro- 
teins were identified. These proteins are indicated to be 
involved in multiple physiological processes, including 
active oxygen removal, defense signal transduction and 
metabolism regulation. PAL and SAMS were up-regu- 
lated at both the protein and mRNA levels. It is clear 
that PAL plays a key role in induced disease resistance, 
and SAMS can catalyzes the production of many resis- 
tance-related material. They can be strongest candidates 
for further functional research. 

Materials and methods 

Plant materials, growth conditions and treatments 

Seeds of soybean (Glycine max [L].Merr.cv. Nannong 
99-10) were sown in 15 cm plastic pots that contained a 
sterile soil mixture (top soil:sand:vermiculite, 3:2:1). The 
seedlings were grown in a greenhouse maintained at 30 
± 5°C, 70 ± 10% RH, and a 14 h/10 h (day/night) photo- 
period with supplementary metal-halide illumination. To 
prevent the interference of other insects, all potted 
plants were covered with gauze. The potted plants were 
watered once daily with the same volume of water. The 
developmental stages of soybean plants were defined as 
described by Fehr et al. [69]. 

Treatments of cotton worm (obtained from a colony 
maintained on an artificial diet) feeding and untreated 
(control) plants were applied to soybean plants at the 
V4 (vegetative stage with four nodes) stage of growth. 
Defoliation (about 20-30%) was induced by placing three 
third-instar cotton worm larvae on the first trifoliate leaf 
for 48 h, after which the larvae were removed. The 
leaves of control and treated plants at four sampling 
times (0 h, 6 h, 12 h and 24 h) after removal of the lar- 
vae were excised for identification of induced resistance, 
protein extraction and RNA extraction. There were 
three biological replications. 

Identification of induced resistance 

To evaluate soybean resistance induced by cotton worm, 
a dual-choice test and force-feeding experiment were 
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conducted. In the dual-choice test, three to five leaves 
from each pair of experimental plants were arranged 
opposite each other at the margins of a porcelain dish 
(30 cm x 20 cm). Ten cotton worm larvae, previously 
starved for 24 h, were released into the middle area of 
the porcelain dish and allowed to feed for 12 h. Leaf 
area was measured by a LICOR-3000 area meter (LI- 
COR, Lincoln, NE). The consumed area in the control 
(C) and treatment (T) was used to calculate the feeding- 
preference index (PI), where PI = 2 T/(T + C) [70]. The 
PI values ranged from 0 to 2, with PI = 1 indicating no 
feeding preference for either control or treatment leaves, 
PI > 1 indicating preference for treatment leaves, and PI 
< 1 indicating preference for control leaves. In the 
force-feeding experiment, control and treatment leaves 
were placed in 15-cm-diameter petri dishes, then five 
cotton worm larvae were placed in the dish. The relative 
growth rate (RGR) was calculated as RGR = (Wl - W2)/ 
Wl, where Wl and W2 are the weights of five cotton 
worm larvae before placement in the dish and after 
feeding for 24 h, respectively. Statistical analysis of the 
data was carried out using Microsoft Excel 2000. Each 
experiment was repeated thrice and the mean values 
and SD were calculated. Single factor ANOVA was car- 
ried out using SAS. 

Total protein extraction 

Total protein extraction for 2-DE was modified from the 
literature [39]. Frozen soybean leaves (0.2 g) were 
ground with liquid N 2 in ice-cold acetone (2 ml contain- 
ing 10% w/v trichloroacetic acid 20 m DTT and 100 
mM PMSF) after incubation at-20°C for 1 h, the pro- 
teins were precipitated by centrifugation at 19,000 x g 
for 15 min at 4°C. The pellet was washed with 3 ml ice- 
cold acetone containing 20 mM DTT five times by cen- 
trifugation at 19,000 x g for 15 min at 4°C until no 
chlorophyll was present in the acetone. The protein pel- 
let was dried in vacuo and resuspended in lysis buffer (7 
M urea, 2 M thiourea, 4% w/v CHAPS, 65 mM DTT, 
0.2% v/v carrier ampholyte [pH 3.0-10.0] and 1 mM 
PMSF) by sonication. The insoluble debris was removed 
by ultracentrifugation at 45,000 x g for 15 min. The 
protein concentration was determined by the Bradford 
method with BSA as the standard [71]. The protein 
samples were stored at -80°C until use. 

Two-dimensional electrophoresis 

The 2-DE experiment was performed essentially accord- 
ing to the Bio-Rad instruction (USA). To proceed with 
isoelectric focusing (IEF), Protein (250 ug) was resus- 
pended in lysis buffer containing 0.2% (v/v) carrier 
ampholyte (pH 3.0-10.0), and 0.01% bromophenol blue 
and rehydrated with each 17 cm, pH 3-10 nonlinear 
immobilized pH gradient (IPG) trip for 13 h. IEF was 



carried out with a Protein IEF Cell (Bio-Rad) at 20°C in 
a stepwise manner: 250 V for 0.5 h, 500 V for 0.5 h, 
1000 V for 1 h, 2000 V for 1 h, 8000 V for 5 h and 
8000 V for 55,000 V h. After IEF, the IPGs were equili- 
brated by shaking at room temperature in equilibration 
buffer (6 M urea, 2% [w/v] SDS, 30% [v/v] glycerol, 50 
mM Tris-HCl [pH 8.8]) containing 1% (w/v) DTT for 
20 min followed by equilibration in buffer containing 
3% (w/v) iodoacetamide for another 20 min. The equili- 
brated IPGs were transferred to a 12% polyacrylamide 
gel. Electrophoresis was performed in Tris/glycine/SDS 
buffer on a Multiphor system (Amersham Pharmacia 
Biotech) according to the manufacturer's recommenda- 
tions. For calibration, low-molecular weight marker pro- 
teins (Amersham Biosciences) were applied on the gel 
via a small piece of filter paper. After electrophoresis, 
the 2-DE gels were stained with silver [72]. The gels 
were fixed in 30% (v/v) ethanol, 10% (v/v) acetic acid for 
30 min, and sensitized for 30 min in 30% (v/v) ethanol, 
0.2% (w/v) sodium thiosulfate and 6.8% (w/v) sodium 
acetate, and then rinsed thrice in water, for 10 mimutes 
for each wash. Impregnated the gels with 0.25% (w/v) 
silver nitrate for 20 min, and then rinsed them twice for 
1 min for each wash in water. 

Transferred the gels in developer solution containing 
2.5% (w/v) sodium carbonate and 0.02% (v/v) formalde- 
hyde until the adequate degree of staining, then trans- 
ferred them to 1.46% (w/v) ethylenediaminetetraacetic 
acid disodium salt stop solution for 30 minutes, last 
washed the gels twice in water for 10 min for each 
wash. 

Gel image comparison and data analysis 

Spot detection, spot measurement, background subtrac- 
tion and spot matching were performed specifically after 
silver- staining of the gels using PDQuest software. Fol- 
lowing automatic spot detection, gel images were care- 
fully edited. Before spot matching, one of the gel images 
was selected as the reference gel. The amount of a pro- 
tein spot was expressed as the volume of that spot 
which was defined as the sum of the intensities of all 
pixels that made up the spot. In order to correct for 
variability owing to silver-staining and to reflect the 
quantitative variation in intensity of protein spots, the 
spot volumes were normalized as a percentage of the 
total volume of all spots present in the gel. The result- 
ing data from image analysis were transferred to 
PDQuest software for querying protein spots that 
showed quantitative or qualitative variations. Only those 
with significant changes (quantitative changes more 
than two-fold in abundance) were used for further 
study. Statistical analysis of the data was carried out 
using Microsoft Excel 2000. The standard deviation 
(SD) was calculated from three spots in different gels. 
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Table 2 Primer sequences used for qRT-PCR 



Spot no 


Protein name 


Forward primer (5'— »3') 


Reverse primer (5'— »3') 


1 


serine/threonine kinase related protein 


1 — rCCCC A A C ! — FTC A A C AT A C A A 

G 1 GGCCAAGC 1 1 CAACA 1 AGAA 


CCl 1 1 1 1 / ' II | r~ r — r A Tl — rA '1 TATAT A 

GGG 1 1 1 1 G 1 1 1 GG 1 A 1 G 1 A 1 1 GA 1 AGA 


2/9 


ribulose-1 ,5-bisphosphate carboxylase/oxygenase large subunit 


CC CC ! — rATTTATTTr A 1 — YC ATfAT 

CGCGG 1 A 1 1 1 A 1 1 1 GAG 1 CAGGA 1 


TV — YCC ( — rC AC A C C A C C C A~T A 

1 G 1 GGG 1 GAGAGGAGGGA 1 A 


3 


Vegetative storage Protein a 


C C A C A A A/"A C C C C 1 — T A A /" A C A 

GGACAAACAGGCCG 1 AACAGA 


t — vr — veer — rc r — vt i i i 1 1 i a~tc a a 

G 1 G 1 GGG 1 GG 1 G 1 1 1 1 G 1 A 1 GAA 


A 


cytoso ic phosphoglycerate kinase 


Trrrr a ata a a tc c cc a *t a 

1 1 GCGAAGAAA 1 GGCGAC 1 A 


c a a a c cc a ct — rrc A C A~TCC A 

GAAAGGGAGG 1 IGAGAIGGA 


c; 


A I r I I LI I die Ufcr Ld iUIJU I II L 


U^UL^ 1 1 nLunLn 1 1 1 O 


C CTTC C A CC AT A CTTC A C CTTC A 

1 1 UUnvJLn 1 r\\J 1 1 unuu 1 1 \Jr\ 


5 


malate dehydrogenase 


IGCCCI 1 1 1 IGCIGAIGCI 


TGCCCCAAGCCCAGAAC 


6 


phenylalanine ammonialyase 


TGCTCAAGGTTGTTGATAGGGAGTA 


TCCACAAGCACTTGCCTTAGC 


7 


S-adenosylmethionine synthetase 


GAGACATGCACCAAGACCAACA 


ATTTCGCGGCATGTGTCA 


8 


predicted protein 


AACAAGTGCAAGTTCGAATCAATG 


ACGACGGCGTAGAGATCGAT 


10 


ascorbate peroxidase 


CGCTCCTCTAATGCTCCGTTT 


AATCGGCGTAGCTCAAAATAGG 


10 


20S proteasome alpha subunit A 


1 ICI 1CGI IGCGLI 1 1 1 ICC 


CUCGAUCAI 1 1 1 IGCCCIAA 




tubulin 


GGAGTTCACAGAGGCAGAG 


GACTTACGCATCACATAGCA 



Gel images used for statistical analysis were obtained 
from three independent biological repeats. 

Protein identification and database searches 

Silver-stained protein spots were excised, washed with 
50% v/v acetonitrile in 0.1 M NH4HCO3 and dried in a 
vacuum centrifuge [73]. Gel fragments were reduced in 
20 ul of 10 mM DTT with 0.1 M NH4HCO3 for 45 min 
at 55°C. After cooling, the DTT solution was replaced 
with 55 mM iodoacetamide in 0.1 M NH 4 HC0 3 . After 
washing, the dried gel pieces were rehydrated in 10 ul 
digestion buffer containing 25 mM NH4HCO3 and 12.5 
ng/ml trypsin (Promega, Madison, WI, USA) and incu- 
bated at 37°C overnight. Tryptic-digested peptides were 
extracted. Samples were analyzed by MALDI-TOF-MS 
and tandem TOF/TOF MS on a time-of-flight Autolex 
III mass spectrometer (Bruker Daltonics, Bremen, Ger- 
many). Peptide mass maps were acquired in positive 
reflection mode, averaging 400 laser shots per MALDI- 
TOF spectrum (resolution was 15,000-20,000). The Bru- 
ker calibration mixtures were used to calibrate the spec- 
trum to a mass tolerance within 0.1 Da. 

The peak list file containing the m/z ratios of precur- 
sor ions and MS/MS fragmented ions was used in 
searches of the web-based search engine Mascot http:// 
www.matrixscience.com against the most recent data- 
base in the National Center for Biotechnology Informa- 
tion (NCBI). The initial search used green plants as the 
taxon. The other parameters for the searches were: 
enzyme trypsin; one or two missed cleavages; variable 
modifications of carbamidomethyl (Cys), oxidation 
(Met), and pyro-Glu (Nterminal Glu); peptide tolerance 
of 0.1-0.5 Da; MS/MS tolerance of 0.2 Da; peptide 
charge of 1+ 2+; and monoisotopic. Only significant 
hits, as defined by the MASCOT probability analysis (P 
< 0.05), were accepted. 



Real-time quantitative RT-PCR 

Total RNA was isolated from frozen soybean leaves by 
the TRIzol method with a total RNA Kit (Tiangen, 
China). Reverse transcription reactions were performed 
with an oligo (dT) 18 primer and ReverTra Ace Moloney 
murine leukemia virus reverse transcriptase (Toyobo, 
Japan) according to the manufacturer's instructions. All 
primers were designed based on the sequences obtained 
using BLASTN in the NCBI database or soybean gen- 
ome sequences database http://www.phytozome.net. 
Real-time qRT-PCR was performed on an ABI 7500 real 
time PCR system (Applied Biosystems, USA) using the 
SYBR Green Realtime Master Mix (Toyobo, Japan). 
Tubulin (GenBank accession no. AY907703), a constitu- 
tively expressed gene with approximately equal PCR effi- 
ciency in all samples, was used as the reference gene. 
Data were analyzed with SDS 2.0 software (Applied Bio- 
systems). The primer sequences used for real-time qRT- 
PCR are listed in Table 2. 
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